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Abstract 

We consider birth-and-death processes of objects (animals) denned in 7L d having 
unit death rates and random birth rates. For animals with uniformly bounded 
diameter we establish conditions on the rate distribution under which the following 
holds for almost all realizations of the birth rates: (i) the process is ergodic with 
at worst power-law time mixing; (ii) the unique invariant measure has exponential 
decay of (spatial) correlations; (hi) there exists a perfect-simulation algorithm for 
the invariant measure. The results are obtained by first dominating the process by 
a backwards oriented percolation model, and then using a multiscale analysis due 
to Klein to establish conditions for the absence of percolation. 

KEY WORDS: birth-and-death processes; random environment; backwards ori- 
ented percolation; multiscale analysis; random point processes; random loss net- 
works. 
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1 Introduction 

A number of reasons explain the standing interest in the study of birth-and-death pro- 
cesses. First, they are important probabilistic constructions in their own right, with 
proven potential to generate and test new mathematical approaches. Second, they offer 
the correct framework to formalize and study several important statistical applications. 
Among the more recent ones, we mention point processes (Baddeley and van Lieshout, 
1995; Strauss, 1995; Baddeley, Kendall and van Lieshout, 1996) and loss networks (Kelly, 
1991, Ferrari and Garcia, 1998). Third, they have become useful tools for other math- 
ematical endeavours. Examples of this are Kendall's (1998) perfect simulation scheme 
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and its offsprings (see M0ller, 2001 for a review) and the study of hard-core statistical 
mechanical systems presented in Fernandez, Ferrari and Garcia (1998, 2001). 

While the non-random version of these processes has a well developed theory (see, for 
instance, M0ller and Waagepetersen, 2004), little is found in the literature for processes 
with random rates. Yet, in many situations a process with spatially homogeneous ran- 
domness is a more realistic model than a spatially homogeneous deterministic one. In 
general, the behavior of processes with "strong" disorder is expected to be qualitatively 
different from that of the non-disordered counterpart. In this paper, however, we tackle 
the complementary issue. We determine sufficient conditions for the disorder to be "weak" 
in the sense that it leads to processes exhibiting ergodic properties not far from those of 
the non-random versions. 

We set our processes on a lattice that for convenience is taken to be Z d . The objects 
being born are called animals. They are supported on finite subsets of Z d , which, for 
technical reasons, are assumed to be uniformly bounded in size. Each animal 7 has an 
associated birth rate w 3 ^) parametrized by a random variable J defined on a certain 
probability space. Each realization of J is a random environment. In addition, there is 
an incompatibility relation between animals which is also assumed to be finite-range. The 
birth-and-death process is defined as follows: for each fixed environment, each animal 
7 attempts to appear at rate w 3 ^) but to succeed, it must past a test, in general of 
stochastic nature, involving incompatible animals present at the moment of the attempt. 
Once born, animals disappear at a unit rate. Our main result, Theorem 12.151 establishes 
two conditions on the disorder [hypotheses ()2.16|) and ()2.18j) below] under which the 
process is ergodic and space-time mixing for almost all random environments. As in other 
models evolving in the presence of frozen-disorder (Klein, 1995; Gielis and Maes, 1996), 
while space mixing remains exponential, time relaxation is only proven to be faster than 
any power, with a bound of the form exp[— mln 9 (l + £)], for some m, q > 0. 

Two conditions are imposed on the disorder. On the one hand, in (12.16)) the average 
local rate of attempted births is asked to have some finite logarithmic moment. This is 
a relatively mild condition which is independent of the incompatibility relation. Its role 
[see Lemma 15.921 and formulas (|5.81|) - (|5.82|) ] is to ensure that there is a not-too-small 
probability that locally no animal is found to be alive throughout a sufficiently thin time 
slice, even in the absence of incompatibility restrictions. The second disorder condition 
()2.18|) requires a weighted birth-rate of incompatible animals to have a sufficiently small 
mean. Weights are defined by a certain size function ^(7) which in principle can be 
chosen in any way leading to the validity of the condition. This is a situation very much 
in the spirit of cluster-expansion formalisms; see, for instance, Dobrushin (1996). Usually 
—but not optimally — the size function is chosen as the number of sites in the support 
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of the animal. This second disorder condition guarantees that large fluctuations in rate 
values are sufficiently sparse to have a negligible effect at sufficiently large scales. 

In a subsequent corollary ( Corollary I2.28JI . we determine a convenient sufficient cri- 
terium for both disorder hypotheses. We show that, barring very exceptional settings, 
there exists an e (which we do not try to optimize) such that if 



E 



sup^|iJ(0)| w J {6)] < e (1.1) 
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the process is almost-surely ergodic and mixing. Here |i?(7)| is the cardinality of the 
halo, which is the region around 7 that determines incompatibility with other animals 
[see (E2U)]. 

Our analysis is based on a construction of birth-and-death processes described in 
Fernandez, Ferrari and Garcia (1998, 2001). In this approach, the process is constructed 
by resorting first to a free process obtained by turning off incompatibilities and allowing 
every attempted birth to succeed. This free process corresponds to a marked Poisson 
process, where the marks are the lifetime of the animal and a random variable to be used 
in the compatibility test. It is useful to visualize each realization of this free process 
as the collection of cylinders determined by the animals alive during a certain lifespan. 
These cylinders inherit the incompatibility relation of the animlas forming their sections. 
Of particular importance for the construction of the interacting process are those incom- 
patible cylinders that are alive when a given cylinder is born. They are called ancestors 
of the latter. The interacting process can be constructed if each cylinder has only a fi- 
nite number of generations of ancestors. In this case, one can descend down genealogical 
trees deciding, by means of succesive incompatibility tests, which cylinders of the free 
process remain for the interacting one. Formally, the relation "being ancestor of" defines 
a (backwards) oriented percolation model in the space of cylinder configurations. The 
construction scheme succeeds if this model does not exhibit infinite percolation. 

A sufficient condition for this lack of percolation is, in turns, obtained by resorting to 
a dominant multitype branching process whose subcriticality allows the construction to 
work. The subcriticality condition is, precisely, 

SU P sT^j E S(9)w(0) < 1. (1.2) 

^ 6 incomp. with 7 

The left-hand side of this inequality is the mean number of branches of the dominating 
branching process, and acts as a driving parameter: Time and space rates of convergence 
and mixing rates can be explicitly obtained in terms of it. 
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This approach can not be directly applied for random birthrates because, except in 
trivial cases, condition (jl.2|) is violated with probability one. An additional argument is 
needed showing that these violations are so sparse that the process is basically driven by 
the behavior in the overwhelmingly present "good" (regular) regions. As often in the study 
of disordered systems, this additional argument takes the form of a multiscale analysis. 
In this paper we adapt a time-tested multiscale argument whose origin can be traced to 
von Dreifus' (1987) dissertation. The argument was later adapted by Campanino and 
Klein (1991), and Campanino, Klein and Perez (1991) to the study of (d + l)-dimensional 
systems with <i-dimensional disorder, which is our setting here. In the present work we 
follow the more general version due to Klein (1994) which has also been the basis of 
Gielis and Maes' (1996) study of spin-flip dynamics. Our proofs are similar to those of 
these references, and the informed reader will recognize many points in common, such as 
conditions ()2.17|) . (|5.43jl and (|5.44|) below. But there is a number of small but frequent 
adaptations needed to accomodate "blob" rather than edge percolation. For the sake of 
clarity and completeness we have prefered to do a self-contained exposition of the proof, 
rather than a catalogue of differences with respect to preceding references. The paper 
size is comparable in both cases. We have, nevertheless, kept most of the notation, and 
of the architecture of the proof, adopted by Klein (1994) and Gielis and Maes (1996) for 
the benefit of readers familiar with these papers. 

The initial ingredient of the multiscale scheme is a sequence of boxes of increasing 
size. The linear sizes of the boxes are the scales. Sites are then classified as regular or 
singular at a given scale according to whether there is an appropriate decay of the inte- 
grated connectivity function to the boundary of a box around the site. For the procedure 
to succeed there must exist a choice of scales that noticeably improves the probability of 
a site to be regular at succesive scales. This requires the singular regions to be sparsely 
distributed within a sea of regular sites; the distribution becoming more and more di- 
luted as the observation scale grows. In the present setting of space-time (backwards) 
percolation, boxes are space-time parallelepipeds with a cubic spatial section. The space 
scales grow as a power law, but the time scale needs to grow faster. This is a conse- 
quence of the (i-dimensional nature of the disorder. While regular and singular regions 
alternate in spatial directions, they are frozen in the time direction. Their effects are 
limited only by the deaths and lack of births of the concerned animals, which are pro- 
cesses of a more correlated nature than the (spatial) disorder distribution. Hence, the 
connectivity function decays much more slowly in time than in space. The time scale 
must satisfy two complementary requirements. First, it must grow fast enough to ensure 
a faster-than-power decay in the percolation probability from bottom to top inside a box. 
But, second, the growth must be suficiently slow to guarantee that time-like connections 
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be mostly established inside a box, rather than by creeping through a sequence of spa- 
tially consecutive boxes. Mathematically, the first requirement is used in the proof of 
Theorem 15.291 to arrive to formula ([5.38)1 . while the second one appears explicitly in the 
proofs of Sub lemma 15.631 and Lemma 15.921 The compromise is achieved by chosen a func- 
tion T(L) so the bottom-top and left-right percolation probabilities inside a box of size 
L x T(L) decrease at a comparable rate, namely exponentially with L. This corresponds 
to a stretched-exponential dependence [(see (j5.114|) - (j5. 116)1 . 

In statistical mechanical terms, our results are, somehow, the low-temperature coun- 
terpart of those of Gielis and Maes (1996). While these authors study single-spin-flipping 
evolutions leading to high-temperature invariant measures, here we focus on gases of "def- 
fects" characteristic of low-temperature invariant measures. The comparison is not fully 
valid, however, because we are able to study only families of uniformly bounded deflects. 



2 Definitions, examples and results 
2.1 Basic definitions 

Animals We consider the lattice Z d , d > 1, (or, in general, the vertices of a graph with 
uniformly bounded coordination number) endowed with some norm. For concreteness we 
adopt = max{|xj| : i — 1, 2, . . . , d}, x G Z d . An animal model in Z d is defined by a 
countable family G of objects, the animals, for which there exists a map 

G — V<?) (21) 

such that (i) V{^y) is a finite set — whose elements are called the vertices or sites of 7 — and 
(ii) there is a finite number of animals asociated to each fixed set of vertices: V" 1 ^^)) 
is finite for every 7 G G. Animal configurations are elements £ of N G . 

Geometry is introduced through the map V in the obvious manner: The diameter 
of an animal 7 is the diameter of its set V(7), the distance between two animals is the 
distance between the corresponding sets of vertices, etc. To abbreviate we shall often 
denote x G 7, for x G Z d , rather than x G ^(7). For a region AcZ d (of finite or infinite 
cardinality), we denote Ga the set formed by all animals 7 with ^(7) C A (animals in 
A). The corresponding configuration will get a subscript A: £ G N Ga . The omission of 
a subscript A indicates A = Z d . Animals are also characterized by a size function. In 
general terms, a size function is any funcion 5* : G — > [l,oo) that can be used in the 
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convergence condition ()2.18|) below. In the examples of next section, this function is just 
the number of sites in V(j). 

In this paper we suppose that there exists a number i\ such that 

diam (7) < £ x for all 7 G G . (2.2) 



Animal interactions We introduce an interaction junction M( ■ | • ) : G x N G — ► [0, 1], 
such that M(7|£) is the probability that an attempted birth of 7 G G actually occurs 
when the current configuration of animals is £. If this function takes only the values or 
1 we refer the interaction as deterministic. This happens, for instance, for conditions such 
as volume- or perimeter-exclusion. The function M determines the (binary) interaction 
matrix 

lnt( 7 ,0) := l{sup|M(7|e)-M( 7 |e + ^)| 7^0} (2.3) 

where 1{^4} is the indicator function of the set A, Sg is the configuration in N G having 
only 9 present, and £ + 5$ is the configuration obtained by adding the animal 9 to the 
configuration £. We say that 7 is incompatible with 9, and denote 7 7^ 9, iff Int(7, 9) = 1; 
otherwise we say that 7 is compatible with 9. A family of animals T C G is a compatible 
family [herd in Dobrushin's (1996) nomenclature] if its elements are pairwise compatible. 
We assume that there exists a finite £2 G R such that 

dist( 7 , 9) > £ 2 => Int( 7 , 9) = . (2.4) 



The random environment Each animal 7 has an associated birth rate w 3 ^) which 
depends on a variable J, the random environment, belonging to a certain probability space 
(J7", P). We assume that the disorder measure P is such that 

T compatible =^ {w 3 ^) : 7 G T} independent. (2.5) 



Birth-and-death processes in random environment For each fixed environment J 
and each A C Z d we consider the interacting birth-and-death processes formally defined 
by the generator: 

AlF(v) = E{ M W ^(7)^ + ^7)-^)] + v(l)[F(v-^)-F(r])}}, (2.6) 

7€G A 

where F is a real continuous function on N Ga . 
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In words, ()2.6|) says that when the current configuration of animals is rj, each animal 
7 appears at rate M(^\rf) w 3 ^) and disappears at rate 1 (if it is present). This is an 
interacting birth-and-death process of animals with quenched disorder. The factor w 3 ^) 
acts as the rate of an internal Poissonian clock marking possible birth instants. The 
actual birth takes place only if a further test, determined by M{^\ ■ ), is passed. From the 
definition (|2.3|) of the interaction matrix we see that if Int(7,6 l ) takes the value 1 (resp. 
0) the presence of 9 may have (resp. does not have) an influence on the birth rate of 7. 

Unit death-rates are no loss of generality. In general, our w 3 ^) stands for the ratio 
"birth rate / death rate" of 7. In fact, a similar treatment is possible for interacting-birth 
interacting-death processes. These are processes with a generator as in ()2.fij) but with 
a factor Md ea th(7|?7) w death(7) multiplying the last summand. Our theory can be easily 
adapted if the death rates wreath (7) are bounded below by a strictly positive number, 
uniformly in 7 and 77. 

In general, there may be forbidden cases for which the matrix M takes value 0. An 
extreme example of this is provided by deterministic interactions (M = or 1), such 
as those defining fixed-routing loss networks and some statistical mechanical models (see 
examples in Section I2.2j) . In these cases, the configurations resulting from forbidden 
cases — that is those of the form £ + S e with M(0\£) = — are not acceptable as initial 
configurations for the interacting process (thus, they will not be generated by it). The 
remaining configurations will be referred to as acceptable configurations. 

Further notation The support of a function F on N G is the set 

Supp (F) = jac e Z d : 3 7 3 x , 77 e N G such that F(rj) ^ F{q + 5 7 ) } . (2.7) 

In this paper we reserve the symbol P for the probability measure on the random 
environment, while the combination Q J (perhaps with further embellishments) will denote 
probability with respect to processes defined by a fixed environment J. We shall use 
superindices for the time coordinate of space-time animal configurations: rf^ e N Ga x IBL 
For finite sets A, the symbol \A\ will denote the cardinality of A. For animals, the notation 
I7] will mean |V(7)|. We shall use a capital letter to denote a space-time point. As a 
default, we shall use the corresponding lowercase letter to denote its space component, 
and a subscript to identify its time component, ex. X = (x,tx) 6^x1. 

2.2 Examples 

Some of the processes fitting our framework are the following. 
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Point processes 



In these models the animals are subsets of Z d , hence V(7) = 7. A large variety of point 
processes has been introduced in the literature. Our framework applies to their space- 
discretized version. The area-interaction point processes of Baddeley and van Lieshout 
(1995) involve identical animals defined by the translations of a fixed compact set G C 7L d . 
An animal of the form 7 = x + G is called a grain of germ x. Animal configurations are 

labelled by the corresponding germs through the identification G D {x + G : x G A} < > 

A C Z d . With this identification, the interaction function takes the form 



M{x\A) = F (x + G)n(A®G) 



X8) 



for some F : N — > [0, 1], where A © G = [j ye A {v + G}. The process is attractive if F is 
increasing and repulsive otherwise. The former case corresponds to the penetrable sphere 
model introduced by Widom and Rowlison (1970). 

In addition, the model is specified by internal-clock rates w(x). In the non-random 
version they are usually independent of x. For the disordered version, the independence 
hypotheses ()2.5|) leads to site disorder: 

J = {J x :xE Z d } (2.9) 

for independent random variables { J x }. The internal clock rates should be of the form 

w \x) = w({J y : y E x + G}) . (2.10) 

As further examples of point processes within our framework we mention: (i) the 
Strauss process (Strauss, 1995), which does not involve grains and the interaction function 
depends on the number of pairs of points closer than a fixed threshold r, and (ii) the 
perimeter-interaction process where the interaction matrix depends on the overlapping of 
perimeters of grains. The latter is a particular instance of the generalization propossed 
by Baddeley, Kendall and van Lieshout (1996). Our results apply to the corresponding 
lattice versions with site disorder. 



Fixed-routing loss networks 

These models are defined on a graph Q whose links define possible calls or connections. An 
example is Z d with the usual nearest-neighbor links, or with a more general family of links 
(x, y) G Z d x Z d with vertices x and y not necessarily at distance one. An animal 7 is a 
connected finite subgraph possibly subjected to some further restriction, depending on the 
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model (eg. not having a loop, or forming a closed circuit). The interaction matrix is often 
deterministic, preventing calls that would cause a link to be used beyond a predetermined 
maximal capacity. More generally, the interaction matrix embodies some penalization 
scheme for the multiple use of links. 

The non-random rates w(j) usually decrease with the number of links in 7. A general 
type of disorder satisfying the independence condition ()2.5|) is site-link disorder: 

J = {J x :xeV(G)}\J{J x , y :(x,y)eL(g)} (2.11) 

with all the random variables J x , J XjV being independent. By V(Q) and L(Q) we denote, 
respectively, the set of vertices and of links of the graph Q. Condition (|2.5|) is satisfied by 
rates of the form 

W J ( 7 ) = w({J x : x G V( 7 )}, {J x , y : (x,y) G L( 7 )}) . (2.12) 

We emphasize that, due to our hypothesis (|2.2|) . our disordered loss network must involve 
connections not exceeding a radius i\. 



Models motivated by statistical mechanics 

These are models of random geometrical objects coming from the study of spin systems 
in statistical mechanics. The associated birth-and-death processes have the statistical 
mechanical measure as invariant measure, and can be seen both as a convenient tool to 
study such measure and as a feasible simulation algorithm for it. We mention two of these 
models. Due to requirement (|2.2|) . the disordered version considered here correspond to 
a "chopped" variant of the corresponding models. 



Random cluster model This model plays an essential role in the study of the Potts 
model. We refer the reader to Grimmett (1995) for its detailed study as well as for 
references to the original articles. The animals, called clusters, are all finite subgraphs 
of a given graph Q (whose vertices are usually the sites of Z d , but whose links may 
include non-nearest-neighbor pairs). The interaction matrix forbides the appearance of 
two clusters having a common vertix. Furthermore, the process has internal-clock rates 
of the form 

»w - n (r^t) n (r) < 2 - i3 > 

for appropiate parameters J X)V and J x [often denoted p(x, y) and q(x)]. A disordered ver- 
sion, satisfying the independence condition ([2.5)1 . is obtained by turning these parameters 
into independent random variables [site-link disorder (j2.11|) ]. 
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Ising contour model This is the model resulting from mapping into Peierls contours 
the typical spin configurations of one of the low-temperature pure phases of the Ising 
model. Contours are hypersurfaces formed by (d — l)-dimensional unit cubes centered at 
points of Z d and perpendicular to the edges of the dual latice Z d + (|, • ■ ■ , |). These cubes 
are called plaquettes. Two plaquettes are adjacent if they share a (d — 2)-dimensional face. 
A contour is a finite family of plaquettes such that (i) the family can not be partitioned 
into two subfamilies with no adjacent plaquettes, and (ii) every (d — 2)-dimensional face 
is covered by an even number of plaquettes of the family. Geometrically, a contour 
corresponds to a connected closed (hyper)surface. The set of vertices ^(7) is the set of 
centers of the plaquettes forming 7. 

The contour-model equilibrium measure is the invariant measure of a birth-and-death 
process with a deterministic interaction matrix which prevents the appearance of two 
contours with adjacent plaquettes. The clock rates take the form 

w { n )= e -^^vM J{x \ (2.14) 

where the parameter (3 is interpreted as inverse temperature and each J(x) is a coupling 
constant associated to the plaquette centered at x. In the disordered version, these 
constants are independent random variables [site disorder (|2.1U|) ]. 



2.3 Main result 

Theorem 2.15 Consider an animal model as defined above, in particular satisfying X2.2]) . 
\2.4\) and Assume the disorder satisfies 



N := E 



ln a (l + sup^w J (7) 



< 00 



73a; 



for some 



a > 2dr 1 + 1/1 + - + 



1 



d 2d 



(2.16) 



(2.17) 



Then there exists an e = e(d, a, K) > 0, with a monotonically decreasing dependence 
K — > e(d, a, N), such that if there exists a function S : G — ► [1, 00) with 



(ii) 



E 



7 GG 5(7) 



< e, 



(2.18) 
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then there exist constants m > 0, q > such that for ¥-almost all configurations J the 
following is true. There exist constants C 3 . > 0, x G Z d , such that: 

1. Existence and uniqueness. For each AcZ d there exists a unique time-invariant 
process Q A on N Ga x R with generator \2.b)) . The process has a unique invariant 
measure /x A . 

2. Time convergence. For each A C Hj d and each acceptable configuration £ . € N Ga ; 

there exists a unique process Q J ( • | £ ) on N Ga x R + defined by the generator \2. b]) 
and the initial configuration £ . Furthermore, the process converges at a superpoly- 
nomial rate to the measure /x A on N Ga . Explicitly, the following inequality holds for 
any function f on N Ga : 



sup 



-=-J 



< 



C J (Supp (/)) exp -m ln"(l + t) (2.19) 



with C J (Supp(/)) = SxeSupp (/) ■ The supremum is taken over all acceptable 
configurations in N Ga . 

3. Space convergence. As A — > Z d ; /x A converges weakly to fi 3 := yU^ d at an exponen- 
tial rate. More precisely, for any function f depending only on animals contained 
in a finite set A: 



n 3 f - vif 



< 



xGSupp(/) 



Cl exp 



-m dist (x, A c ) 



(2.20) 



4- Exponential mixing. For any functions f and g depending on animals contained 
in an arbitrary set A C Z d : 



A(fg) -AfA.9 



< 



\g\\oo /2 ^lCy 



exp 



-m\\x-y\\/2 . (2.21) 



xGSupp (/), 
J/GSupp (g) 



5. Perfect simulation. The invariant measure /i J ( or any of the /x A j can be perfectly 
simulated using the ancestors algorithm introduced by Fernandez, Ferrari and Garcia 
(2002). 
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These properties are similar to those of the deterministic process obtained by a fixed 
translation-invariant asignment of J. The major difference is the possible slowing-down 
of the time convergence due to the subexponential time dependence in the right-hand side 

of (irroji . 

We remark that while the size function 5* involved in hypothesis (|2.18|) is completely 
arbitrary, it is clear that it conveys some idea of the "mass" or "might" (Dobrushin's, 
1996, terminology) of animals, in the sense that the larger Si^y) the larger the set of 
animals incompatible with 7. In fact, the natural choice for the examples presented in 
Section 12.21 is the size of the relevant set of vertices of each animal (area of the grain, 
length of the call, perimeter of the contour, etc.). But, of course, "natural" is not a 
synonym of "optimal" . 

In most cases of interest, condition ()2.18|) is morally more limitant than ()2.16|) . To 
formalize this fact, let us denote 

T J = su P ;>> J (7) (2.22) 

x 

and 

* J = -g^)E^)» J W, (2-23) 

and let us introduce the halo of an animal, namely the region of 7L d around it beyond 
which no incompatible is possible: 

H(i) := p| \\V C Z d : 7(0) C W 9 ~ 7}' . (2.24) 

This halo is a finite set because of assumption (J2.4j) . In the examples of Section 12.21 
if (7) = 7(7). In addition we denote 

\Hh)\ \Hh)\ 
u x = inf U2 = sup^-^ (2.25) 

7 eG: 7 ^0 S(7) 7eG 5(7) 

and define 

S J = sup^|#(0)| w 3 (9) . (2.26) 

x e^x 

In terms of these quantities we have the bounds 



which lead to the following corollary of our main result 



T J < S J and ^ J < — S J (2.27) 

Ul 
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Corollary 2.28 Consider an animal model as defined in Section 

u\ > and w 2 < oo . 



such that 



(2.29) 



Then, there exists e > such that z/E(H J ) < e there exist constants m, q > such that 
properties 1 to 5 of Theorem \2.1® hold for ¥ -almost all disorder configurations J. 

The conditions ()2.29|) are trivially satisfied by all the examples in Section 12.21 for 
which, in fact, ui = u 2 = 1. In our bounded-diameter setting, they are also satisfied if 
the function size is chosen translation invariant. 

Proof. We first observe that for any a > 1 there exists C a such that ln a (l + x) < C a x 
for all x > 0. Therefore, 



E 



ln a (l + T J : 



< C n E 



(T J ) 



Pick an a as in (J2.17)) and take 



mm 



{77- > —e{d,a, 1)) , 



(2.30) 



(2.31) 



where e(d, a, K) is the function given in Theorem 12.151 From ()2.30|) and the leftmost 
inequality in (|2.27j) we obtain 



E 



ln a (l + T J ) 



< C n e < 1 . 



(2.32) 



From the rightmost inequality in (J2.27j) and the monotonicity of the function K 
e(d, a, K) we obtain 



E[^ J ] < — e < e(d,a,l) < e(d,a,$) . 



(2.33) 



Thus, both hypothesis ()2.16|) and ()2.18|) of Theorem 12.151 are satisfied. □ 



3 The basic construction and the key lemma 

The proof of Theorem ()2.15j) relies on two ingredients: a graphical construction of birth- 
and-death processes propossed by Fernandez, Ferrari and Garcia (1998, 2001), and a 
multiscale argument following Klein (1994). In this section we review the graphical con- 
struction and state the key result obtained from Klein's argument. 
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3.1 Graphical construction 

The construction is performed for each fixed environment J. 

The free process We associate to each animal 7 an independent marked Poisson pro- 
cess iV 7 with rate w 3 ^). We call ^(7), 7 G G, the ordered time-events of iV 7 with 
the convention that To (7) < < ^(7). For each occurrence time ^(7) of the pro- 
cess iV 7 we choose independent marks 5* (7) exponentially distributed with mean 1 and 
Zii^f) uniformly distributed in [0,1]. The free animal process is the process in which 
at each Poisson time-event ^(7) (a copy of) the animal 7 appears and lasts ^(7) 
time units. It is convenient to identify each marked point (7, ^(7), 5*. (7), ^(7)) with 
(7 x \T k ('j),T k ('y) + Sk(j)] , Zfc(7)), the cylinder with basis 7, birth time ^(7) lifetime 
Skd) and mark Z^j). This identification turns the free process into a measure on the 
space of cylinders. Denoting C = (7, t, s, z), we use the notation 

Basis (C) = 7, Birth (C) = t, Life (C) = [t, t + s], Mark (C) = z 

The marks will be used later on to define the (interacting) birth-and-death process with 
generator (|2.6jl . Let us denote Q J the probability measure on Z d x R corresponding to 
the free process (it always exist, being a countable product of marked Poisson processes). 

Backwards oriented percolation We first extend our definition of incompatibility 
to cylinders in the natural way: Two cylinders C and C are incompatible if they have 
incompatible bases and they are simultaneously alive at some instant of time. That is, 
C r/, C if and only if Basis (C) ^ Basis (C) and Life (C) n Life (C) ^ 0; otherwise 

Let us now fix a family of cylinders C_ (for instance, obtained as a realization of the 
free process of the previous paragraph). We define the ancestors of a cylinder CeCas 
the set 

Af = {C e C : C </> C and Birth (C) < Birth (C)| . (3.1) 
Recursively for n > 1, the nth generation of ancestors of C (in C_) is 

A^ = \C" : C" G Af for some C" G A^_ x } . (3.2) 

[In fact, A^ = A„ (C), but we shall omit this dependence except if it is crucially needed.] 
Likewise, for an arbitrary space-time point (x, t) G 7L d x K we define its set of ancestors 
as the set of cylinders that contain it 

A^' t] = |(7 G C_ : Basis (C) 9 x, Life (C) 9 *} , (3.3) 
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and, recursively, 



Ai*>*> = \C" : C" E Af ' for some C E A^l } • (3.4) 

Definition 3.5 We say that there exists backwards oriented percolation in C_ if there 
exists a space-time point (x, t) such that An ^ for all n, that is, there exists a point 
with infinitely many generations of ancestors. 

The clan of a space-time point (x,t), resp. of a cylinder C, is the union of the corre- 
sponding ancestors: 

A (..*) = (J A^) , A c =[jA n c . (3.6) 



n>l n>l 



Other quantities that will be used later are the time-length, space-diameter and space- 
size of the clan of a point (x, t): 

TL (A (x '* } ) = t - infjs : Life (C) 3 s, for some C G A ( ^} (3.7) 
SD(A^) = diam( |J Basis (C)) (3.8) 

SS(A (x '* ) ) = j [J Basis (C) . (3.9) 



The interacting birth-and-death process If 

Q J ({no backwards percolation}) = 1 , (3.10) 

the stationary process with generator ()2.6|) can be constructed by "cleaning" the free 
process defined above. For completeness we present a summary of this construction here. 
The reader is referred to Fernandez, Ferrari and Garcia (2001) for details. The idea is to 
start from first ancestors ( "Eves" ) and classify cylinders into kept or erased according to 
the test determined by the interaction function. Cylinders that are born in presence of a 
kept ancestor and that fail the test are erased, all the others are kept. 
Explicitly, let C_ be a cylinder configuration in the set 

C = {C_ without backwards percolation} . (3.11) 

Since all clans in C_ are finite, each cylinder has a well defined, finite number of ancestors. 
Therefore the configuration can be decompossed in the form C_ = U n > (2 n , where C_ n : = 
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{C G C : A^ 7 7^ 0, A^ +1 = 0}. The sets if and .D of kept and erased cylinders are defined 
inductively as follows. Starting with Kq = C_ and D = (cylinders without ancestors 
are kept, as they do not need to pass any test), we define, recursively, 

K n = {CeC n \ U-- \D t U JQ : Mark (C) < M(C\ U^ 1 JQ)} 

(3.12) 

R n = C n \[^uu^ui I )] 

where M(C|^') = M(Basis (C)\ {Basis (C) : C" G Birth (C) G Life(C')}) . We 
denote the set of kept cylinders as K_(C) = U n K n and the set of erased cylinders as 
D{C) = U n D n . Clearly 

(CI) C_ is the disjoint union of K_(C_) and D(C_), and 

(C2) The event {C G K(C)} is measurable with respect to the sigma field generated by 
A c , in fact K_(A C ) = K_{C_) H A* 7 . In words, it is sufficient to know the (finite) clan 
of C to know whether C is kept or erased. 

The stationary animal process Q J is defined by the sections of the kept cylinders: 

rf'\l,C) = l{ Basis (C) = 7, Life(C) 3 t} . (3.13) 

CdK{C) 

If Q J (C) = 1, the process Q J is Markovian and has generator (|2.fij) : that is, 

|E J /(^(C)) = E J Aif{rf{C)) (3.14) 

This fact is proven in Fernandez, Ferrari and Garcia (2001) for the homogeneous case. 
This proof extends to the inhomogeneous case in an obvious manner. Let us denote fi 3 
the distribution of any t-section rj J,t , which, by construction, is independent of t. We shall 
determine the properties of /i J by studying the law of rj 3,0 , the stationary (interacting) 
birth-and-death process at time zero. 

As observed in (C2) above, the presence or absence of contours intersecting a region 
V C 1j d at time t depends only on the clans of the cylinders alive at time t whose bases 
intersect V, that is, on 

A y >* := jC G A c : Basis (C) D V ^ 0, Life (C) 3 t} . (3.15) 
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In particular the function ?? J '*(7, ■ ) defined by (|3.13|) is in fact a (deterministic) function 
only of A 7,i [= A^ 7 ^*]. More precisely, if C_ and C' are two cylinder configurations such 
that A 7 '*(C) = A 7 '*(C'), then r/ J '*(7,C9 = ^ J '*(7, C). We code this fact as the identity 
(slightly abusive from the notational point of view): 

77^(7, C) = 7^(7, A 7 ' 4 ). (3.16) 



3.2 The key lemma 

Let us call a sequence of cylinders Ci, C2, . . . , C n an open path if C2 76 Ci, . . . , C n 7^ C n _i 
and Birth (Cj+i) < Birth (Cj) for all i. Given a cylinder configuration (7 and two space- 
time points X = (x, t) and Y = (y, s), with s < t, we say that X and K are connected (in 
the configuration (7) if there exists an open path C±, C2, . . . ,C n such that x G Basis (Ci), 
t G Life(Ci), ?/ G Basis (C,,,), s G Life(C„). [Equivalently, K is in (the interior of a 
cylinder belonging to) the clan of ancestors of X.} The existence of such a connection 
defines an event denoted X — > K. 

For a given realization of the environment, the connectivity function is defined by 

G 3 (X, Y) = Q J {X — > F} . (3.17) 

Many of the properties stated in the main theorem are a direct consequence of the 
following result 

Lemma 3.18 Under the hypotheses of Theorem \2.1fft [namely, \2.2\) . \2.J$ , \2. 5}) and 
\2.11\j ], but replacing \2. 1 6}) by an inequality 



E 



ln a (l + sup^w J ( 7 )) < N, (3.19) 



79x 



there exist m > and qo(a, d) > 1 such that for allq G (1, go) there is a value e(d, a, m, q, N) > 
so that the validity of \2.1ty) implies that for every x G 7L d 

G J ({x,t),{y,s)^ < Cl expj-mmax ||x - y\ , ln 9 (l + \t - s\) } (3.20) 

for all y G Z d and t, s G K, where the constants C^. = C%(£i + h) are finite for F-almost 
every environment J . 

The fact that the inequality ()3.19|) is used implies that the function K — > e( • • • , N) 
can be chosen to be decreasing, as stated in Theorem 12.151 
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We will prove this lemma by performing a multiscale analysis similar to the one used 
by Klein (1994) in his work on extintion of contact process in a random environment. 
That will be done in Section We first discuss how Theorem 12.151 follows from the 
bound (J3.2(Jj) for the connectivity function. 



4 How the theorem follows from the key lemma 
4.1 Relation with percolation properties 

The following theorem relates properties of the measure fi 3 with properties of the perco- 
lation model. 

Theorem 4.1 Assume that for a given J there is no backwards oriented percolation with 
Q 3 -probability one. Then, 

1. Uniqueness. The measure fi 3 is the unique invariant measure for the process r) 3,t . 

2. Time convergence. For any function f with finite support, 

4/-E J (/(^) | i h 



sup 



< 2 



to J {TL(A^°)) > bt\ +e-( 1 - fc ) t E J (sS(A^°)))j (4.2) 



xeSupp (/) 



for any b G (0, 1). 



3. Space convergence. As A — > Z d ; /i^ converges weakly to [i 3 . More precisely, if 
Supp (/) C A ; then 



V 3 f - t*lf 



< 2 



xeSupp (/) 

4- Mixing. For f and g with finite support, 



Q J {SD(A( a,0 >) > dist({x},A c )}. (4.3) 



Vi(fg)- m Va9 



< 4 



9\ 



oo y oo 



[Q J {SD (A( 2 '°)) > ||x - ?/|| /2} + Q J {SD (Afo'°>) > \\x - y\\ /2} 



a;GSupp(/), 
yeSupp (g) 



(4.4) 
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Proof. Items 1, 2 and 3 follow from displays (4.6), (4.7) and (4.9) of Theorem 4.1 of 
Fernandez, Ferrari and Garcia (2001). The analogous of item 4 is stated in that theorem 
with the extra assumption that there exists a time h such that there is no space-time 
percolation in (0, h). We provide here a proof without this hypothesis. 

We consider functions / and g such that Supp (/) = {x} and Supp (g) = {y}, the 
general case follows by telescoping. Let's fix a partition {T, T'} of Ga- Below we choose 
r formed by animals "closer to x" . Let A and B_ be Q J -independent realizations of {C G 
C : Basis (C) G T} and A[ and B[ independent realizations of {C G C : Basis (C) G P}. 
Then A U A, A U £>', B U A', Bug have the same law as C A and A U & is independent 
of B_ U A!. Let A^ X '^(C) be the random variable defined in (|3.6|) and 

X(A,A!) = /(^(-iUA')) = /(^(..A'^UUi'))). 

Analogously 

^01,4') = ff^V.iUi')) = g(t| J,0 (-,A M (AuA'))). 
With these definitions, we obtain 



/4(/<?)-/4//4<? = E x(A^')^(A^)-^(A5')r(5,A') 



(4.5) 



where E corresponds to a four-fold product of the measure 
to the bound 



This expression leads us 



< 2 



with 



A = [X(A,A!) = X(A,g) and Y{A,A) = Y(B, A) | . 



(4.6) 



(4.7) 



We now choose V as the set of animals intersecting {z G Z d : \\z — x\\ < \\z — y\\} and 
T' as its complement. Then the event A is verified whenever the bases of the cylinders 
in both A^ X '°\A U A!) and A^(A U g) are contained in V and those of A^°\A U A) 
and AW(BUi') are contained in P*. The complement of the intersection of these four 
events yields 

P(i c ) < 2Q J |SD(A (:E ' 0) ) > \\x-y\\/2^ + 2 Q J {sD (A (i/ ' 0) ) > ||x - y|| /2} . □ 



For comparison purposes, let us present an alternative mixing bound. 
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Proposition 4.8 Assume that for a given J there is no backwards oriented percolation 
with Q J -probability one and consider functions f and g with finite support. Then, 



< 2 



\9\ 



Q J {SD(A (a; ' 0) ) + SD(A fe ' 0) ) > ||ar-y||} 



(4.9) 



xeSupp (/), 

j/GSupp (g) 



where Q J is the free process obtained doubling the birth rates of Q J . 

The bound ()4.9|) corresponds to standard high-temperature results in statistical me- 
chanics (see, for instance, the main result in Bricmont and Kupiainen, 1996). Its proof 
relies on the very popular technique of "duplication of variables". In contrast, (J4.4j) is 
proven by "tetra-plication of variables". In our general setting, however, we are able to 
exploit better our first bound ([4.4)1 . 

Proof. As above, it is enough to assume Supp (/) = {x} and Supp (g) = {y}. Let C_ and 
C' be Q J -independent realizations of Ca- We denote 

X(C) = f(v 3 '°(-,C)) and Y(C) = g(r] J % ,£)) . 

The duplication-of-variables identity is: 

1 



Ma(/<7) - Ma./>a<? = 2 E 
where E corresponds to the measure 



X{C) - X(C')j (Y{C) - Y(C') 
x O j . Let us now consider the event 



B = {Basis (A^GUC')) f] Basis (a^ 0) (C U C\ 



(4.10) 



(4.11) 



and the transformation T : (C,C!) — > (TC,TC!), that interchanges the (x, 0)-ancestors 
in C and C: 



TC 
TC' 



C\A^°\C)j UA {X '°\C') 
C!\A^(C')) UA^(C) 



(4.12) 



Conditioned to B being true, the distribution of (C_,C') coincides with that of (TC_,TC!) 
(the processes inside and outside each realization of Basis ^A^ X '°\C U (7') J are indepen- 
dent). Furthermore, the event B is T-invariant. Hence, 



E 



F(C,C!) 



E 



F(TC, TC!) 



(4.13) 
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for each local function F on Ga x Ga. But, in the presence of £>, the function F involved 
in (j4.10j) is odd under this transformation. We conclude that 



E 



(x(C) - X(C'fj (Y(C) - K(C7')) 1{B} 







which, by (|4.1()jl . implies 

v>l(fg) ~Kf K9 



< 2 



I^HooP(0 C ). 



(4.14) 



(4.15) 



The proof follows from the observation that 

F(B C ) < Q J {SD (A (x ' 0) ) + SD (A fe ' 0) ) > ||a: - y\ 



□ 



4.2 Bounds for the size of the clan 

The key inequality (|3.20|) leads to the following bounds for the probabilities of the time- 
length and space diameter and size of the clan of a space-time point. 

Proposition 4.16 Under the hypotheses of Theorem \2.1bl there exist m,m>0 and con- 
stants C x , and C^., x G Z^, which are finite for F-almost every environment, such 
that 



J {TL (A (x ' 0) ) > t} < C J X expj-m ln 9 (l +T)} 
> J {SD(A( X '°)) > l\ < Cl exp{-mL} 



E J {sS(A (a: ' 0) )} < C 3 X . 
Proof. As max (a, b) > (a + b)/2, inequality ()3.20|) leads to 

G J ((x,t),(y,s)} < C x J exp{-^[||x-y||+ln 9 (l + |t- S |)]} 
Inequality (J4.17)) is then a straightforward consequence of the fact that 
Qj | TL(A (x-,o ))>r | < J2G J ((x,0)^(y,-T) 



(4.17) 
(4.18) 
(4.19) 

(4.20) 

(4.21) 
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while ()4.18|) follows from 




(4.22) 



To obtain (|4.19|) we use (|4.18|) to bound the right-hand side of the inequality 

Q J |SS (A (x ' 0) ) > l} < Q J {SD (A (x ' 0) ) > (4.23) 

and sum over L. □ 

4.3 Proof of Theorem 12^51 

The bound ()3.20|) provided by the key Lemma 13.181 implies, by Borel-Cantelli, the ab- 
sence of backwards oriented percolation. Thus, we can apply Theorem 14.11 for almost-all 
realizations of the disorder. 

1. Existence and uniqueness. Under the hypothesis of no percolation in any re- 
gion A, the process ()3.13|) defines a stationary process rf in A by considering those 
cylinders in C with basis contained in A. Uniqueness follows from the next item. 

2. Time convergence. We use ()4.17|) and ()4.19|) to bound the right hand side of ()4.2|) 
by 

211/lloc (ptew{-™W{l + bt)} + die-W) 

zeSupp (/) 

< 2H/IU C 3 x exp{-m\tf(l + bt)} } (4.24) 

zeSupp (/) 

for some constant C% < oo. 

3. Space convergence. Inequality ()2.20j) follows immediately from ()4.3j) and (|4.18|) . 

4. Exponential mixing. The bound (|2.21|) is just the combination of (|4.4|) and (|4.18|) . 

5. Perfect simulation. The construction proposed in Fernandez, Ferrari and Garcia 
(2001) works, for almost all disorder realizations, due to the absence of backwards 
oriented percolation. 
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5 Proof of Lemma 13.18 



The lemma is proven by adapting Klein's (1994) multiscale analysis. The analysis involves 
several steps that will be studied separately. 

5.1 General scheme and notation 

The scheme is based on a sequence of linear sizes, called scales, defining boxes of increasing 
size. Sites are classified on regular or singular according to the decay of connectivity 
functions on a surrounding box. The definitions must be tunned up so that the probability 
for a site to be regular increase sufficiently fast with the scale. The presentation of this 
scheme is organized as follows: 

(i) In Section I5~2l we recopilate a number of inequalities needed to prove next inductive 
step. In particular, the Hammersley-Lieb-Simon inequality is the key to pass from 
local to global decay of the connectivity function. The uniform bound on the size 
of the animals is needed to ensure its validity. 

(ii) In Section 15.31 we define regularity and determine "good enough" probabilities for 
sites to be regular at each scale. Regularity at all scales with such probabilities 
guarantee, by Borel Cantelli, almost sure regularity. 

(iii) Section E31 shows the crucial inductive step: A "good enough" probability of being 
regular at a given scale implies a "good enough" probability of being regular at the 
next scale. This is the heart of the argument. 

(iv) The last step of the proof is the determination of conditions so the origin has a 
"good enough" probability of being regular at some initial scale. This is done in 
Section 15.51 

We introduce some notation. We fix £q = £± + £2, the maximal size of the animals plus 
the radius of incompatibility - to avoid trivialities we assume £q > d + 1. For L > and 
iGZ" 1 , we denote 




(5.1) 



and define the £ -boundary of A[x; L] as the set 



d io A[x;L] = A[x; L + 5] \ A[x; L), 



(5.2) 
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where 5 > 1 is chosen so the following is true: For any sequence of animals 71, 72, • • • , 7n 
with 7j 7i+i, 2 = 1, 2, . . . , n — 1 connecting a point inside A[x; L] with a point outside 
A[x; L+8], at least two animals in the sequence are contained in dg A[x; L\. This condition 
is necessary to satisfy the Hammersley-Lieb-Simon (see next). A possible choice is 5 — 

Up \ _ 3(lp-2) 

<H*oj - 2(3=1)" ■ 

For X = (x, t) G Z d x R, L > and T > we set 

S L)T (X) = A[x;L] x[t-T,t}. (5.3) 

The vertical, horizontal and complete boundaries of the box B LjT (X) are defined respec- 
tively as: 

d v B L , T (X) = d to A[x;L] x [t-T,t], (5.4) 
d H B L>T {X) = A[x;L + S] x {t-T}, (5.5) 
dB L:T (X) = dvB LjT (X)Ud H B L:T (X). (5.6) 

Note that, as we are considering backwards oriented percolation, the face A[x; L + 5] x {t} 
is excluded from the boundary. 

For any integrable function H : Z d x R — > R we denote 

]T i?(Z) = ^ H(z,t-T)+ f H (^ s ) ds - ( 5 - 7 ) 

ZedB L , T {X) zeA[x;L+<5] zefy,A[x;L] *~ T 

We introduce the notion of connection within a box. Given a cylinder configuration 
C_, its restriction to a box B = A x I , with AcZ d and / a real interval, is the family £7 B 
of cylinders obtained by "restricting" to / those cylinders of C_ with bases inside A. That 
is, each C = (7, t, s,z) G C_ with 7 G Ga defines a cylinder Cb = (7, t/, sj, z) G C s with 
tj = max(t, inf I) and sj = min(t + s, sup /) — ij. Given two space-time points, X = (x, t) 
and Y = (y,s), X,Y G B with s < t, the event X — > F is formed by all cylinder 
configurations C_ such that the configuration C_ B exhibits an open path connecting X 
with Y. If this event is true, we say that X and Y are connected in B. The event defines 
the connectivity function in the region B: 

G J B (X,Y)=Q 3 {X^ B Y}. (5.8) 

As usual, we shall omit the subscript if B = Z d x R. We write 

g 3 b l+sMX) (x, d) = Cm ( x > z ) • ( 5 - 9 ) 
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5.2 A toolbox: Inequalities 



The first two inequalities we need have been basically proven in the literature. The minor 
adaptations needed for our setting do not justify a detailed exposition of their proofs. We 
content ourselves with providing appropriate references and indications. 

We consider the natural partial order in the space of cylinder configurations: C_<C^ 
if C[ contains all the cylinders in C_. Events are said increasing, resp. decreasing, if their 
characteristic functions are nondecreasing, resp. nonincreasing with respect to this partial 
order. 

Proposition 5.10 (Harris-Fortuin-Kasteleyn-Ginibre inequality) If A and B are 

both increasing or both decreasing events, 



In the context of continuous-time percolation, an inequality of this sort was first proven 
by Bezuidenhout and Grimmett (1991). They did so by showing that the process is the 
weak limit of discrete-time processes that satisfy the corresponding inequality. Their 
approach provides at the same time a proof of the van den Berg-Kesten inequality, but 
has the inconvenience of imposing an additional topological condition on events, namely 
that their boundary — in some suitable metric topology — have measure zero. While such 
requirement is satisfied for the events of interest to us, a different approach, based on the 
martingale convergence of expectations of the relevant events, yields the result as stated 
above. A proof along this line is provided in the book by Meester and Roy (1996) (see 
their Theorem 2.2) for the percolation of spheres of random radii on a continuous space 
(Poisson Boolean model). The proof is easily adaptable to our setting of percolating 
cylinders. 

The second inequality refers to increasing events happening in a disjoint manner. For 
brevity, we state it only for the type of events needed in the sequel. Let us consider a 
box B = A x I, with AcZ d and I a closed finite interval in R, and space-time regions 
Bi, . . . ,B n C B. We denote {B x — >b B 2 } o • • • o {B n _i — > b B n } the event of having 
n — 1 open paths, respectively connecting in B some point in Bi to some point in B i+ i, 
1 < i < n — 1, such that no two paths share the same cylinder. 

Proposition 5.12 (particular case of the van den Berg-Kesten inequality) ForB, 
Bi, . . . ,B n as above, 



Q 3 (AnB) > Q 3 (A)Q 3 (B) . 



(5.11) 




'n-l 



B fin}} < Q J { 



}---Q J { 



B n -\ >b B n 



(5.13) 
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This inequality is a consequence of the more general inequality proven by van den Berg 
(1996) for increasing events of marked Poisson processes. The only subtlety is that van 
den Berg's result requires the events to depend on Poisson clocks ringing within a bounded 
region of R d+1 . The connectivity events in (J5.13)) do not seem to fit in this framework 
because they may be determined also by cylinders born in an arbitrarily remote past. To 
obtain (|5.13|) we must, therefore, apply van den Berg's result to the marked process in 
B obtained by adding, independently, an invariant initial distribution of cylinders. This 
corresponds to a further, independent, spatial Poisson marked process on A x {inf 1} 
with rates {if -1 (7) : 7 G Ga}- Less direct proofs are also possible either by adapting the 
Bezuidenhout and Grimmett (1991) approach mentioned above, or the proof of Theorem 
2.3 of Meester and Roy (1996). 

Our last inequality is obtained in Klein (1994) as a corollary of van den Berg's in- 
equality. In our oriented setting, we can present a totally different argument. 

Proposition 5.14 (Hammersley-Lieb-Simon inequality) Let L,T G M+, 5 as de- 
fined below and W C Z d x R. Then, for every X = (x,t x ) G W and Y = (y,t Y ) G 
W \ Bl + s,t{X) , with ty < tx, we have 

G 3 W (X,Y) < G k +S , T (X)(X,Z)G 3 W {Z,Y Z ) (5.15) 

z edB L:T (X)nw 

with Y z = {y,mm(t Y ,t z ))- 

Proof. The enclosing character of Bl + $^t{X) and the choice of 5, implies that ev- 
ery connection from X to Y includes a connection within Bl + ^^{X) to a point Z G 
dBL,T{X) fl W. It is here where the assumption of bounded animal sizes is critically 
needed. This point Z is subsequently joined by an open path to a final cylinder whose 
section contains y. If ty <t z , this open path determines the event Z -^> w Y (recall that 
only connections that are backwards in time are considered). If on the contrary tz < ty, 
then this final cylinder must have been born before t z and stood alive at least up to ty. In 
particular, it contains the point (y,tz). Both situations are summed up in the inequality 

G 3 W (X,Y) < Yl ® 3 { X ^b l+sM x)Z; Z ^ w Yz) . (5.16) 

zedB L:T (x)nw 

(The sum in the right-hand side involves, in fact, a time integral. Its justification requires, 
therefore, a limit of conectivities involving time-discretizations of dyB^^iX). These 
connectivities are continuous functions, so the limit converges to the integral. We omit 
the details.) 
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The desired inequality ()5.15|) follows from ()5.16|) and the fact that the events X — > b l+s t (x) 
Z and Z — > w Y z are independent. Indeed, they are determined, respectively, by which 
cylinders are alive for t < tz and which are alive for t > tz- Such events are independent 
due of the exponential character of cylinder lives. Perhaps the simpler way to see this is by 
resorting to an alternative construction of the cylinder process, where for each animal 7 
the birth and death events are generated independently with respective exponential rates 
u> J (7) and 1. Cylinders are born at a birth event and live up to the next death event, 
neglecting intermediate events of the wrong type. With this construction, connectivity 
events before and after a given time t z are determined by different, hence independent, 
birth and death events. □ 

Corollary 5.17 Let L,T 6 R + , 5 as defined below UTty and W C Z d x R. Then, for 
every X = (x,tx) and Y = (y,ty) in W with ty < tx, we have that, 

G 3 W (X,Y) < g b l+s , t (x)( X ' Z i) G b l+s , t (z 1 )( z i, Z 2 ) ■■■ 

ZiedB L , T (x)nw z 2 edB L , T (Zi)nw 

^ Gbl+s^Zn^^n-i, Z n ) G 3 w {Z n , Y Zn ) 

Z N edB L , T {Z N -i)C\W 

(5.18) 

for every 

N < integer part of max ( — — ^ , ■ X Y - ) . (5.19) 

V L + T J 

Proof. This is just an iteration of (|5.15jl . The number of times such iteration can be 
performed is at least equal to the right-hand side in (j5.19j) . □ 



5.3 Regularity and "good enough" probabilities 

In this section we introduce the main notions defining the multiscale approach. 

From now on we consider boxes as in (|5.3jl where the temporal height T is an increasing 
function of L. In Section T5.4I the function T(L) will be eventually chosen as a stretched 
exponential, but the following results do not depend on such a particular choice. To 
simplify the notation we characterize the boxes by its spacial lenght L and denote 

B L (X)=B L+mL) (X). (5.20) 

Following Klein (1994) we separate Z d in regular and singular regions for a fixed 
realization of the environment J. 
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Definition 5.21 Let m > and L > 1. A site x G Z d is said to be (m, L)— regular if 

G j Bl{x>o) ((x,0) , d) < e~^ L+5 \ (5.22) 

Otherwise x is called (m, L)— singular. A set A G Z d is said (m, L)— regular if every 
x G A is (m, L)— regular; otherwise it is (m, L)— singular. 

Regularity will be used in conjunction with the Hammersley-Lieb-Simon inequality 
through the following crucial result. 

Lemma 5.23 Let A be a (m, L)-regular region, W D A x R. Then, for every X = (x, tx) 
and Y = (y, ty) in W , 

G J W (X,Y) < exp{-m(L + 5)N A (X,Y)} (5.24) 



with 

\ x -y\\ \tx-t Y \ 



' dist(x, A c ) 

A^a(A^, Y) = integer part of min <{ — - — ^— — , max 

L + o 



L + 5 ' T 



(5.25) 



Proof. This follows from Corollary 15.171 The value of N\(X,Y) satisfies the constraint 
()5.19j) and guarantees that all the intermediate sites Zi are in A and, hence, they are reg- 
ular. We then bound the right-hand side of (|5.18jl starting from the right: Gw(Zn,Yz n ) 
is bounded by one, and each of the preceding sums by exp{— m (L + 5)}. □ 

We also formalize the notion of scale 

Definition 5.26 A scaled sequence is a triple (Lo,a,T) where Lo,a > 1 and T : 
(0, oo) — > (0, oo) is a function that grows faster than any power. Each such triple defines 
an increasing sequence of sizes L k+ i = L% for k = 0, 1, 2, . . .. The length L k is the k-th 
scale of the sequence. 

Finally, we associate "good probabilities" to scales. 

Definition 5.27 A scaled sequence (Lo,a,T) has moo-good-enough probabilities if 

there exists p > ad such that 

pjx is (rricojLk) — regular j > 1 — — p (5.28) 
for all k = 0, 1, 2, . . for all x G Z d . 
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We end this section with the proof of the "easy part" of the multiscale argument. 

Theorem 5.29 If a scaled sequence (L ,a,T) has m^-good-enough probabilities, then, 
for any m G (0,moo) there exist constants {C^{m) : x G Z d } with P{C^(m) < oo : x G 
Z d } = 1, such that 

G 3 (X,Y) < Clim) exp j -m max \\x - y\\ ,T- l (t x -t Y ) } (5.30) 
for all X, Y G Z d x E with t x > t Y G R. 

Proof. We follow Klein (1994). Take b > 1 (to be determined later) and consider 

S k '■— |a[x; b(L k+ i + 5)] is not a (moo, L k ) — regular region j. (5.31) 

This event is verified if at least one of the sites in A[x; b(L k+ i + 5)} is (m^, L fc )-singular. 
As, by hypothesis, the probability for a given site to be singular is at most L^ p , we obtain 

Pjy is (moo, £&) - singular | 

i/GA[a;;6(L fc+1 +5)] 



< 



[26(L 



fe+i 



6)] d ^ {AbLtY 



LI 



< 



(46) 



- a k (p—ad) 
J 



As p > ad, this bound shows that the probabilities P{5fc} are summable in k. Therefore, 
by Borel-Cantelli, with probability 1 there exists k\ = ki(x,£o,b,J) < oo, such that 
A[x; b(L k+1 + 5)] is (m^, L k )- regular region for all k > k±. 
Now fix X, and classify the sites Y into regions 



n 



Y : max 



x-y\\ , T-\t x -t Y ) < b(L kl +5) 



(5.32) 



and 



|r: b(L k + 5) < max \\x - y\\ ,T-\t x -t Y ) < b(L k+1 + 5)} (5.33) 



for k > k\. 

Let us first consider Y G lZ k for some k > k\. In this case we have ||x — y\\ < 
b(L k+ i + 5) and so y G A[x;b(L k+ i + 5)} which is a (moo, L fc )-regular region. It follows, 
from Lemma (for W = 7L d x R and A = A[x; b(L k+1 + 5)]), that 



G\X,Y) < e 



-m ca (L k +S)N 



(5.34) 
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with 



v • t , t ■ )H L k+l + 5) 
JM = integer part of mm < , max 

L k + o 



\z-y\\ t x -ty 



L k + 5' T(L k 



(5.35) 



A first bound of (|5.34j) comes from the observation that, as Y e lZ k implies \\x — y\\ < 
b(L k+1 + 5), 

'\x - y\\ 



N > integer part of 



L k + 5 



Therefore, 



G 3 {X,Y) < exp\- moo (L k + 6) 



If - y\\ 

L k + 6 



(5.36) 



(we bounded the integer part of a number by the number minus one). 

This bound can be improved in cases where the temporal part dominates in the sense 
that T _1 (tx — ty) > \\x — y\\. In this case we use that b(L k + 5) < T~ l (tx — ty) for 
Y G TZf., so that, as T grows faster than any power, 

tx-ty > T(b(L k + 6)) > b(L k+1 + 5) > \\x - y\\ 



T{L k ) ~ T(L k ) L k + S ~ L k + 5 

if hi, and hence k, is chosen large enough. Hence 

» r s • + , b(L k+1 + 5) fT'^tx-ty) 
JM > integer part of — - — — - — > integer part of — — (5.37) 



L k + 5 



L k + 6 



where the last inequality comes from the definition of lZ k . From this and (|5.34p we obtain 
that 



G 3 (X,Y) < exp { - moo (L k + 6) 



T-\t x -t Y ) 
L k + 5 



(5.38) 



whenever T 1 (tx — ty) > \\x — y\\. 

Inequalities ()5.36|) and ()5.38j) can be combined in the expresion 

G 3 (X, Y) < expj-moo [max(||x - y\\ , T~\t x - ty)) - {L k + *)]}■ 



(5.39) 



The additive correction (L k + 5) can be turned into a factor (1 — 6 x ) because (L k + 5) < 
(1/6) max{||x — y\\ , T _1 (tx — ty)}, by definition of lZ k . In addition we choose 6 = n J n °^ m 
so that m = m^l — r). In this way ()5.39j) yield 



G J (X,Y) < exp j-m max \\x-y\\ , T _1 (t x -t F ) J 



(5.40) 
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uniformly in k > k\. 

Finally we consider Y G 1Z. They satisfy max{||a; — y\\ , T _1 (ix — ty)} < b (L^ + 5), 
so we can write 

G 3 (X,Y) < G 3 (X,Y) exp|m6(L fcl expj-mmax ||ac -y\\ , T'^tx -t Y ) } 

< C 3 (£ ,m) expj-mmax ||a; - y|| , T _1 (tx- - t y ) J (5.41) 
where C 3 (£ ,m) = e mb( - Lk i +5 \ The desired bound ijOJI) follows from (EHUJ) and (fCTjl . □ 

5.4 The change of scale 
5.4.1 The change-of-scale theorem 

This section contains the heart of the multiscale argument. Its main result is the following 
theorem, related to Theorem 3.2 of Klein (1994), which establishes sufficient conditions 
for a scaled sequence to have good-enough probabilities. 

Theorem 5.42 Let d > 1 and choose a satisfying \2.1T\ i. Assume that the disorder is 
such that K < oo satisfies \3.19i) for the chosen a. Let a = d + y(P + d and take v and p 
such that 

ad [a + a + 1] 

—T i ;r < v (5.43) 

a [a — d + ad\ 

and 

a{v\a — d + ad] — ad) — ad , 

ad < p < . (5.44) 

a 

Finally, consider m andm^ such thatO < < m . Then, there exists L = L(d, 5, a, K, u, p, mo, rn^) < 
oo such that if for some L > L 

p|x is (m , £o) ~ regularX > 1 — —p , (5.45) 
L J L 

for all x e Z d , then the scaled sequence (L ,a,e L ") has m^- good- enough probabilities. 

Note that condition ()2.17|) guarantees the existence of v satisfying ()5.43j) . which in 
turns implies the existence of p as in ()5.44j) . 

The combination of this theorem and Theorem 15.291 implies that the key Lemma 181 
— and hence all the properties stated in Theorem 14 .11 — follow once f!5.45|) is satisfied. Note 
that the logarithmic time-dependence in the key lemma is a consequence of the choice 

T(L) = exp(L u ) , (5.46) 
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< v < 1, for the time- height of the scaled sequence. 

The proof of Theorem 15.421 is inductive. We show that it is possible to arrange the 
parameters in such a way that, if 

p|x is (m, I) — regular J > 1 — — (5-47) 

for x G Z d , m > and I sufficiently large, then we also have 



p{ac is (M, L) - regular } > 1 - (5.48) 



where L = l a and M is an appropiate function of m. The inductive step has two ingredi- 
ents. First, in S ubsect ion 15 . 4 . 21 we show that the inductive regularity holds in the presence 
of two events A and Ba- This part of the argument is based on geometric considerations 
and the inequalities of Section 15.21 Subsequently, in Subsection 15.4.31 these events are 
proven to hold with high-enough probability. The success of the approach relies on the 
careful definition of the events A and Ba- 

5.4.2 Geometrical estimates: Good events imply good behavior 

Lemma 5.49 Consider some fixed x G Z* d . Let u, a, I be such that 0<u<l,a>l and 
I > 5. Set 

9 := min{a — 1, a(l — u)} (5.50) 

and take mo, 9, m with m > 0, < 6 < 9 and l~ e < m < m . Put L = l a , pick a positive 
integer R and define the event 

A = | J : there exists Xi,x%, . . . , Xr G A[x; L + 5] such that 

R 

A[x; L + 5} \ A[xj] 2(1 + 5) + 1] is a (m, /) — regular region^ . (5.51) 

3=1 

Assuming A is true, take A > 0, b such that < b < ctv and k > max{l, v + 6*0}. Define 
the event 

B A = {J: II K ^ * e- 1 "), (5.52) 

where 

R 

A = U^A[^ ; r]f|A[a;;L + 5]) (5.53) 
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—the union being taken on the points involved in the definition of the event A — and 

Ki{i) = e -( 1 + A )«' J (7) + (i _ e — J W)(1 _ e ~ A ) e - AwJ ^ . (5.54) 

Then, there exist ao = ao(d, 5, a, v, mo, R) > and Iq = lo(d, 5, a, v, k, mo, 0, 8q, R) < oo 
such that if I > Iq the following holds: If A and Ba are true for some A G (0, 1] then the 
site x G Z d is (M, L) - regular with 

M>m-^>±. (5.55) 

This lemma is a direct consequence of the following Sublemas 15.561 and 15.631 where 
we analyze separately the connectivity of the site X = (x, 0) to sites in the vertical and 
horizontal boundary. The regularity stated below (|5.54|) follows immediately by summing 
the corresponding bounds over both boundaries. 

Sublemma 15.561 involves estimations purely in the spatial direction, hence only the 
event A is relevant and the choice of the time scale T(L) actually plays no role (though, 
for concreteness, we stick to the subexponential dependence). Event £>a is required to 
control the time-like percolation studied in Sublemma 15.631 It provides a lower bound for 
the probability of not having long towers of ancestors based around "deffective" sites [see 
(|5.82j) - (j5.83|) below]. The stretched-exponential choice (|5.46|) becomes essential for this 
sublemma which is in fact the hardest estimation of the paper. 

Sublemma 5.56 There exist li = l\{d, 5, a, v, mo, 9, R) < oo anda\ = a\(d, 5, a, v, mo, R) > 
such that for all I > h, all Y G dyB^X) and all J G A 

G j Bl{x) (X,Y) < e~ M ^ (5.57) 

with 

Ml =m-^>±. (5.58) 

Proof. We work with a fixed environment J G A. First we group the deffect-sets 
A[xj]2(l + 5) + 1] into larger cubes which absorb "pockets" totally surrounded by orig- 
inal deffects. This leaves us with a much simpler geometrical situation, where the reg- 
ular region is the complement of a finite family of cubes. More precisely, elementary 
geometrical considerations show that there exists a possibly smaller collection of sites 
Vx, U2, ■■ ■ ,Ur> G A[x; L + 5], with R' < R, n 1 ,n 2 , ...,n R > G {1,2,. . .,R} and n l + n 2 + 
. . . + tiri < R, such that the sets A[yi] n«(2(Z + 5) + 1)], i = 1,2, ... ,R', are at a dis- 
tance strictly larger than one, [L=i M x j> 2(i + 5) + 1] C Uili MVi'i n «(2(Z + 5) + 1)], and 
A' := A[x; L + 8] \ [jf =1 A-tyi] n>i(2(l + 5) + 1)] is a connected (m, Z)-regular region. 
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These problematic regions define cylinders of radius rij(2(Z + 5) + 1) and temporal 
height e L " centered at the points y{. 

Bi = B n .( 2 (i+s)+i),e L " ((Vi, 0)) , i = l,2,...,R 

[recall the notation in (|5.3J) ]. We shall control the connectivity function in the spacial 
direction, on the (m, Z)-regular region 



R' 



B' = 5 L (X)\|JS 4 



(5.59) 



i=l 



[we use the notation resulting from ()5.20|) and ()5.46j) ]. 

If X G B', we denote dB = {X}; otherwise X G By for some i' and we put dB = 
OBii. Similarly, if Y G B' , we denote dB R+1 = {Y}; otherwise Y G Byt for some i" and 
we put dB R+ i = dBiu . Every connection from X to Y can be decomposed into disjoint 
connections among some of the cylinders Bi. Therefore, 

R' 

r=l {i 1 ,i 2 ,...,i r }C{l,2,-,R'} 

o---o{dB ir ^ B ,dB R , +1 }. (5.60) 
As B' C A' x R with A' a (m, Z)-regular region, we can apply (J5.24j) to obtain 

12 t /Da 



\dL v 



exp < — m(l + 5) 



31,32 



1 + 5 



Q J {dB h >b' dB n } < [(R[2(l + 5) + l])~e 
We have denoted 

Dj lt j 2 = min|||xi — a^H : (^i,t) G dBj 1 , (x2,s) G <9£>j 2 for some t,s| 
and used the bounds 



. (5.61) 



max 



\ x ji x j2\\ r s l 



1 + 5 e L " \ ~ 1 + 5 ~l + 5 
^From (|5.6Uj) . (|5.61|) and the van den Berg-Kesten inequality (|5.13j) we have 

R 1 



Gl L(x) (X,Y) < J2 



dL v 



(R[2(l + 5) + l]) d e 



2(R+l) 



r=1 {ii,*2,...,v}C{l,2,...,fl'} 

J n , x\ I ^0,ii + Ai,i 2 + • • • + A r ,_R'+l / . 1 • 

x exp < -m(£ + oj I — ■ I + 5 ( r + ^ 
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This inequality, together with the bounds D 0}il + Di ui2 + . . . + r'+i > L — 2[2(/ + 8) + 

1] [n h +n l2 + ... + n lr ] > L — 2[2(l + 5) + 1]R and ^Zi Z {il ,i 2> ... M c{l,2,...,R' } 1 < (R + 1) ! > 
yield 

G J BLm (X,Y) < (R+l)<[(R\2(l + S) + l}) d e L "]' {R+11 

x exp {_ m( , + s) ( L-W + S) + W _ (R + 1} ) } . (562) 

Therefore 

G J B L( x)(^n<e-^ (L+5) ; 

with 

»w cim c 2 

M x = m 



for some fixed constants Cx, c 2 > depending on <i, 5, a, z/ and i?. This yields the proposed 
bound (15371) because, as m > y^-, 

M( > M 1 = m-£ > -1 

for some ai = ai(d, 5, a, z/, m , -R) > and for I large enough (the meaning of "large 
enough" depends on d, 5, a, z/, m , 6 and R). □ 

Sublemma 5.63 Suppose the events A and in \5.51)) and are true. Pick r such 

that v < t < min{/t — 6*0, au}. Then there exists l 2 = h{d, 6, a, u, mo, 8, k, b, r, R) < 00, 
such that for I > I2 we have 

G 3 BL(x) (X,Y)<exp{-M 2 e^} 

for allY e d H B L (X), with 

M 2 = m- e~V > -L . (5.64) 

Note that it is possible to choose r satisfying the hypothesis because k > v + 8 and 
v < av. 

Proof. Let us fix an environment J G A D £>a- The proof relies on the introduction of 
additional space and time scales: 

6(1) = l K and r(l) = exp(/ r ) . (5.65) 
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The time scale r(Z) is intermediate between T(l) and T(L) = T(l a ) because 1 < r < av. 
The space scale S is used to construct the region A defined in (J5.53j) . For the time being 
it is only required to grow faster than / (/c > 1), but the hypotheses of Lemma 15.921 
below will force k to be strictly smaller than a [see condition (|5.1U2j) ]. making 5 a scale 
intermediate between I and L = l a . The intermediate time scale defines a partition of the 
box Bl(X) into 

T(L) 

N = integer part of (5.66) 

T{1) 

"slices" of time-heigth r(Z): 

S 3 = B L+sm ((x, -(j - l)f (/))) , j = 1, 2, . . . , N . (5.67) 
We note that events in different Sj are independent (see the paragraph preceding Corollary 
We introduce the notation Ba = A x [— T(L), 0] for any A C Z d . Let 

R 

A = ( |J A[x i; 2(/ + 5) + 1]) p| A[x; L + 5] (5.68) 
i=i 

be the set of "deffective" or "irregular" sites. The region A defined in (|5.53j) corresponds 
to a "buffer zone" around, and including, these irregular sites. We decompose the event 
{X — >b l (x) Y} according to whether there exists at least one slice such that no deffective 
point within it is connected to a preceding slice. That is, we write 

[x ^ Bl{x) y) = ({x ^ Bl{x) y } n a) (J ({x -> bux) y) n A c ) (5.69) 

with 

N N 

A = U A i = U{ b a n 5 i — ^ ^+i} c • ( 5 - 7 °) 

The probability of the first event on the right of (J5.69|) is damped by regularity. Indeed, 
the occurrence of {X —*b l (x) ^} H ^ implies the existence of a connection in B A ^ x . L+g ^ A 
of vertical height at least equal to r(l). Therefore, 

{X^ Bl{x) Y}nA c \j{(yi, Sl ) ^ ^a[x;l+s]\a (y^ s ^} . ( 5 - 71 ) 
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the union being taken over all (yx, si), (2/2, S2) £ -^Ab-i+5]\A "with I s i — s 2 |> t(Z). By 
regularity fLemma l5.23|) we get, after simple computations 



({X - fljl(JC) Y}nA 
< exp< — m(Z + 5) 



m 

T{1) 



- 1 



+ 2d ln(2L + 25) + 2 In T(L) } . 



(5.72) 



Therefore, 



({x -+ BlW y}ha) 

< exp< — m(l + 5) 



e - l 



2dln(2P + 25) + 2/ c 



< exp 



V V 



(5.73) 



for / large enough (recall that v < r). 

To prove that under our hypotheses A c is sufficiently improbable, we bound 



Ai D F;C\D_ 



<j-l/2)?(!) , 



(5.74) 



where Fj is the event that there is no connection inside the slice Sj from points in the 
deflective set to regular points outside the buffer zone A: 



F 'J = { dB A n S 3 >SAB X B A[x;L+5]\A R S 3 } 



(5.75) 



and Z) s represents the lack of connection within a strip of height A involving deflective 
sites: 



A x {s} 



Ax[s-A,s] 



A x {s - A} 



(5.76) 



We are assumming that I is sufficiently large so that t(I)/2 > 1, and hence r(l)/2 > A 
for all A G (0,1]. 

Since each Fj involves connections within regular regions, we can apply Lemma [5.231 
to obtain 

Q J {F<} < e~ h{l \ (5.77) 

with 



h(l) = m{6{l)-(3l + 4d + l)}-2\nr{l)-dlnS(l)-dln(2l + 3S+l)-2\nR. (5.78) 
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Hence 

h(l) > mc x l K -c 2 l T > c 3 r-° > c 3 l K ~ eo , (5.79) 

where ci,c± and C3 are positive constants — depending on d,8,K,T,9 and R — and / is 
large enough. We have used that m > l~ , 9 < 9 and r < re — 9. 
On the other hand, 

D s d P| {[£(7) U f 2 ( 7 )] n £3(7)} , (5.80) 

where 

£1(7) = {7 is n °t present at time s — A} , 

£2(7) = {7 is present at time s — A but it does not survive until time s} , 
£3(7) = { there is no birth of 7 in the interval [s — A, s]} . 

We observe that 

Q J ([^(7)U^(7)]n^(7)) = (q J [£i(7)]+Q J [£ 2 (7)])Q J [£ 3 (7)] 

= e — J (7) + (!_ e — J (7))(!_ e -A) 



e -A» J (7) 



= Kid). (5.81) 

This result follows from well-known properties on Poisson processes and from the fact that 
the events in which an animal 7 is activated at time —t and survives until time define an 
inhomogeneous Poisson process of rate w J (7)e _t . Therefore, due to the Q J -independence 
of events pertaining to different animals, 

Q J {D S } > J] K J A (j) > e- 1 " . (5.82) 
7 ec A 

By the FKG- inequality we get, from 1(5771)1 . (I5~77|) and (l5~79l - (l5~HTl) . 

Q J {^} > Q J {i^} Q J {D_jj_i/2)r(0} 

> (l - e" 

> e- 21 " (5.83) 
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for I large enough. Thus, the independence of the events Aj, j — 1, 2, . . . , N leads to the 
bound 

N 

Q J (A C ) = Jjfl-Q'iAj)) < (l-e- 2l Y < exp(-Ne~ 2lb y (5.84) 
i=i 

The condition r < av implies that, for I large enough, N > exp(Z Qiy /2), and as, besides, 
b < au, we obtain from (J5.84)) that 

Q 3 (A C ) < exp(-e r,//4 ) , (5.85) 

for I sufficiently large. 

Finally, inequalities ()5.73j) and (|5.85j) together with the decompostion ()5.69|) yield 
(recall that r < au) 

G j Bl(x) (X,Y) < Q J ({X^ BL{x) Y}nA)+Q J (A c ) 

( _i t i t ~i r i av ~) 

< exp< — (m — e - ^") > + exp< — e~ > 

< exp{ — (m — e~^)e~\ , (5.86) 



with m — e r ^ 3 > 1/L e , for / larger than a certain threshold that depends on d, 5, a, u, 
m , 9, k ,b, t and R. □ 

5.4.3 Probabilistic estimates: Good events have high probability 

In this section we show that, assuming (|5.47)L the events A and £>a introduced in (|5.51|) 
and ()5.52|) have good enough probability for I sufficiently large. Hypothesis (j3.19|) is used 
to bound the probability of the second event. 

Lemma 5.87 Let m>0,l>5,a>l and p > ad. Put L = l a . Pick a positive integer 
R > apj {p — ad) and define the event A as in \5.51\) . If 



p|x is (m, /) — regular | > 1 — — (5. 



for all x G Z d , then 

PM} > 1 - (5-89) 
assuming I large exceed a certain minimum value that depends on d, a, p and R. 
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Proof. Let us say that two sites £i, 22 G 1i d are Z-nonoverlapping if the distance between 
the boxes A[xi, I + S) and A[x 2 ; I + 5] is strictly larger than one. In this case, the events 
{xi is (m, /) — regular}, i — 1, 2, are independent. We have 



i?+i 



*4 C C [J P| jxj is (m, Z) — singular j, 

x 1 ,...,x R+1 eA[x;L+S\ i=l 



(5.90) 



where the index (/) reminds us that the union is over collections of R + 1 sites which are 
/-nonoverlapping. Hence 



(0 

'} < £ 

x 1 ,...,x R+1 &A[x;L+S] 



P< x is (m, I) — singular 



-iR+i 



IP(R+1) 



LP l(p~ ad ) - 2Lp ' 



(5.91) 



for I large enough. □ 



Lemma 5.92 Let I > 5, a > 1 andp > 0. Put L = l a . For k and b such that < k < b/d 
take a collection of sites {xi : i — 1, 2, . . . , R} and introduce the event £>a as in \5. 5fy) for 
A = e~ lv , with < r] < b — nd. If K < oo for some a > a{p + d)/r), then 



p{#a} > i 



l 



2Lp 



(5.93) 



for I larger than a certain threshold that depends on a, p, k, b, d, rj, K. 

Proof. We start with the following observations: 

(i) For any collection of sites {xi : % = 1, 2, . . . , R} 

R 
i=l 

(ii) As a consequence of (i), 



(5.94) 



B C A C 



U { E ln 

y£A[x;L+8] -f&GA[y;l K +S] 



> 



R 



(5.95) 
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(iii) For any A > 0, 



^a(7) > (l-e- A )e- AwJ(7) > 0. 



(5.96) 



^From ()5.9(ij) and the inequality ln(l — e A ) > ln(A/2), valid for A small enough, we 
obtain 



1 



< [A W J ( 7 )-ln(l-e- A ) 

7GG A [y;iK+5] 



zeA[y;l K +S] 732 



< 4 c WAsup^u> J ( 7 ) -6 G ln(— ) 



(5.97) 



with 9q := sup x |{ 7 3 x}\ < oo. Therefore, for A = e < 77 < 6 — and assuming 
/ large enough in every step we get 



l^ G A[y,l K +S] 



1 



Liri( 7 )J R 



> 



) < p|Asup^ W J ( 7 )>^- + e G ln(-)| 



79a; 



b—nd 



< p{su P j:»-'( 7 )>- ffr _} 



79x 



< p{ln(l + sup^w J ( 7 ) 



> 



79a: 



< — E 



ln a M +supJ^w J ( 7 ) 



(5.98) 



79a; 



where in the last step we have used Chebyshev inequality. Finally, ()5.95|) . f)5.98j) and the 
fact that at] > a(p + d) imply that 



1 o( a+M )K 1 



U> lari-a(p+d) — 2L P 



(5.99) 



for / large enough. □ 



5.4.4 Discussion on optimal choices 

We pause to discuss the possible choices for the parameters involved in Lemmas 15.491 
15.871 and 15.921 The parameter a in Lemma 15.921 involves explicitly the birth rates of the 
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animals. We observe that the smaller the value of a the weaker the condition on the birth 
rates. Thus, the goal is to choose the smallest possible value of this parameter compatible 
with the hypotheses. Let us first list the requirements on the different parameters: 

a > 1 ; (5.100) 
< v < 1 ; (5.101) 

„ b olv 

max{l,i/ + o } < k < - < — ; (5.102) 

< r] < b-Kd; (5.103) 
p > ad- (5.104) 

a > (5.105) 
V 

Inequalities ()5. 102)1 tell us that for choice of k and b to be possible we must have that 
olv > d and hence that [see (j5.5()jl ] 

9 = a(l - v) . (5.106) 

This identity, in combination with the extreme and leftmost inequalities in ([5.102)1 . re- 
spectively implies that 

v > , 5.107 

a — d + ad 

and 

k > v + a(l-v). (5.108) 
Note that (I5.1()7|) yields, in view of (|5.H)lj) . that 

a > d. (5.109) 

On the other hand, if we combine the rightmost inequality in ()5.102|) with f)5. 108|) and 
(|5.103|) we obtain 

i] < u[a — d + ad] — ad < a — d (5.110) 

[the last inequality is due to (jS.lOlj) ]. Finally, we use this last bound, in combination with 
()5.104|) and (|5. 105)1 to obtain a lower bound for the parameter a: 

ad(a + 1) . . , . „ _ _ s 

a > — — ^ with a > d . (5.111) 

a — d 

Small values of a compatible with this restriction are obtained when 

a = d+Vd 2 + d (5.112) 
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and satisfy the inequality 



> 2^(1 + ^1 + - + -). (5.113) 

Reversing the preceding analysis, we see that the best strategy is to take a and a acording 
to (|5.112|) and (|5.113|) and to choose v and p satisfying (|5.43j) and (|5.44|) . This yields the 
smallest possible choices of a that allow choices of rj, b and k respecting the constraints 
(j5.100j) - (j5.105|) . This explains the definition adopted for the parameters in the statement 
of the change-of-scale Theorem 15.421 

This is also a good opportunity to comment on choices for time scales. They are de- 
termined by the demands impossed throughout the proofs of Sublemma 15.631 and Lemma 
15.921 There are two major constraints: First, the scales 5(1) and r(Z) must be truly 
intermediate in the sense that 

I ■< 5(1) -< L = l a and T(l) -< r{l) -< T(L) = T(P) , (5.114) 

where /(/) -< g(t) means that f(l)/g{l) tends to zero as I tends to infinity. Second, the 
exponent h(l) defined in f)5.77j) / (j5.78j) must grow with I. This implies that 

r(l) -< e*® . (5.115) 

The combination of the two preceding displays yields the relations 

T(l) -< f(Z) -< e 5{l) -< e [a (5.116) 

which shows that at best the time scale can grow as a stretched exponential. 

5.4.5 Conclusion of the proof of Theorem 15.421 

We prove (|5.28J) by induction in k. The case k = follows from hypothesis ()5.45j) because 
(mo, L) -regularity gives (moo, L)-regularity for m^ < m . Suppose now ()5.28|) is true 
for some k > 0. Lemmas 15.871 and 15.921 — with the replacements (m, I) — > (mk, Lk) and 
L — > — and the preceding discussion on optimal choices imply that if 

p(x is (m k , L k ) - regular) > 1 - , (5.117) 

then 

P{^lni3 A } > l--=p— (5.118) 
L k+i 
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for I, and hence L k , sufficiently large. On the other hand, Lemma f5 .491 — with (m, I) — > 
(m k , L k ) and (M,L) — > (m fc+1 , L fc+1 ) — implies that if L^ e < m k < m for < 9 < 9 = 
a(l — z/), then 

AhBa C jx is (m fc +i, Lfe+i) - regular j , (5.119) 

with 

(2q 1 / \ 

m k+1 > m k --^ > -g— (5.120) 

for I, and hence sufficiently large. ^From ()5.118|) and ()5.119|) we get that 

Pja; is (m k , L k ) - regular} > 1 ~ Jp ( 5 - 121 ) 

for all natural k, for I sufficiently large. 
To conclude we must check that 

< m k , for k = 1, 2, . . . (5.122) 

^From ()5.120|) we have 

fc-i . fc-i 



m k > m -a J2 7^ = m o ~ a oJ2 (j^) ■ (5.123) 
j=o l j i=o VL o 

Therefore, f)5. 122|) is verified if L is chosen so large that 

a °5Z(7^) <"io-^oo- □ (5.124) 



5.5 A last probabilistic estimate: The choice of initial scale 

To finish the proof of the key Lemma 13.181 we show that for every real L (though we are 
interested in L large) there exists L > L such that ()5.45|) holds. It is only here that 
hypothesis ()2.18|) is needed. Let us suppose that a function S : G — > [1, oo) has been 
selected and consider the resulting quantity \1/ J defined in f!2.23|) . 

Lemma 5.125 Given L,p > 0, for each p 6 [0, 1) there exist m(p),e(p) > and L(p) > 
L such that 

P|^ J > p} < e{p) Pjx is (m(p),L(p)) - regular} > 1 - -— y . (5.126) 
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Proof. By Theorem 5.1 in Fernandez, Ferrari and Garcia (2001), the condition \I/ J < p 
implies that there exist c(p),m(p) > such that 

G b l+s ,t(x)^ Y ) < <P) exp{-m(p) max(||x-</|| ,-t Y )} , (5.127) 

for all Y = (t/jty) with ty < 0. Taking, for instance, m(p) = m(p)/2 we conclude that 
there exists L(p) such that 

jx is (m(p),L) - regular | D < p| (5.128) 

for all L > L(p). Implication (J5.126|) follows by choosing L(p) = max[L, L(p)] and 
e(p) = L(p)-v. □ 

Chebyshev inequality shows that ()2.18|) is a sufficient condition for the validity of 
the left-hand side of f!5.126|) . as long as e < ps(p), for any given choice of p. This 
observation, the preceding Lemma [5.1251 and Theorems 15 .421 and 15.291 constitute the proof 
of the main Lemma lH.181 We observe that, in our proof, 1 = \ and go — ad\a+a+i} w ^ n 
a = d + \J d 2 + d. 
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